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SWIMMERS 



T. CHAMBRION AND A. MUNNIER 

Abstract. In this article, we are interested in studying locomotion strategies 
for a class of shape-changing bodies swimming in a fluid. This class consists of 
swimmers subject to a particular linear dynamics, which includes the two most 
investigated limit models in the literature: swimmers at low and high Reynolds 
numbers. Our first contribution is to prove that although for these two models 
the locomotion is based on very different physical principles, their dynamics 
are similar under symmetry assumptions. Our second contribution is to derive 
for such swimmers a purely geometric criterion allowing to determine wether 
a given sequence of shape-changes can result in locomotion. This criterion can 
be seen as a generalization of Purcell's scallop theorem (stated in [9]) in the 
sense that it deals with a larger class of swimmers and address the complete 
locomotion strategy, extending the usual formulation in which only periodic 
strokes for low Reynolds swimmers are considered. 



1. Introduction 

1.1. About Purcell's theorem. The specificity of swimmers at low Reynolds 
numbers (like microorganisms) is that inertia for both the fluid and the body can 
be neglected in the equations of motion. Consequently, as highlighted by Purcell 
in his seminal article [9 j, the mechanisms they used to swim are quite counter- 
intuitive and can give rise to surprising phenomena, the most famous one being 
undoubtedly illustrated by the so-called scallop theorem. Roughly speaking, this 
theorem states that periodic strokes consisting of reciprocal shape-changes (i.e. a 
sequence of shape-changes invariant under time reversal) cannot result in locomo- 
tion (i.e. does not allow to achieve a net displacement of arbitrary length) in a 
viscous fluid. Considering the prototypical example of the scallop (as sketched on 
the left of Fig. [T]), which can only open and close its shell, and assuming that the 
animal lives in a low Reynolds environment (which it does not), Purcell explains 
that it can't swim because it only has one hinge, and if you have only one degree of 
freedom in configuration space, you are bound to make a reciprocal motion. There 
is nothing else you can do. In addition to the light this result casts on the under- 
standing of the hydrodynamics of swimming microorganisms, it has to be taken 
into account as a serious pitfall for the design of micro-robots, for which engineers' 
interest grows along with the number of applications that have been envisioned for 
them (such as, for instance, drug deliverers in the area of biomedicine). 

The classical assumptions of Purcell's theorem are that the shape-changes have 
to be time periodic and the sequence of shapes (over a stroke), invariant under time 
reversal. Notice that the latter condition does not mean that the shape-changes 
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Figure 1. On the left, Purcell's 2D scallop cannot swim in a vis- 
cous fluid... and neither can it in a potential flow. On the right, 
3D Purcell's corkscrew can swim in a viscous fluid and probably 
also in a potential one. 

have to be strictly time-reversal invariant, with the same forward and backward 
rate, but only that the succession of shapes is the same when viewed forward and 
backward in time. Under these hypotheses, Purcell concludes that the swimmer 
comes back to its initial position after performing a stroke. Going through Purcell's 
article, one will find no proof for this result. However, a huge literature devoted 
to this topic has been produced since then and mathematical proofs can be found, 
for instance, in the article of E. Lauga and T.R. Powers [5] (which contains also an 
impressive list of references and to which we refer for a comprehensive bibliography 
on this topic) and in [2 by DeSimone et al. 

1.2. Beyond Purcell's theorem. Although, as already mentioned, Purcell does 
not provide a rigorous proof of his famous theorem, he explains that the keystone 
of his result relies on that inertia is not taking into account in the modeling of low 
Reynolds swimmers, allowing in particular the Navier- Stokes equations governing 
the fluid flow to be simplified into the steady Stokes equations. Our first main 
contribution in this article will be to prove that more widely, Purcell's theorem 
in its original form still holds true for a class of swimmers subject to a particular 
linear dynamics that will be made precise later on. This class obviously includes 
low Reynolds swimmers but also high Reynolds swimmers extensively studied in the 
literature (see for instance the article [3] of E. Kanso et al. or [1 by T. Chambrion 
and A. Munnier, and references therein). 

Purcell's theorem does not admit any reverse statement allowing to determine 
wether a sequence of shape-changes violating the hypotheses can result in locomo- 
tion. To illustrate this idea, consider Fig. [2] on which is plotted the graphs of the 
functions t G R+ \-> Oj(i) (for j = 1,2), where t stands for the time and each 0j(i) 
(.7 = 1, 2) gives the value of the angle of the scallop's hinge, as sketched on the left 
of Fig. [I] None of these sequences of shape-changes is neither periodic nor time 
reversal invariant. However, anybody familiar enough with Purcell's result would 
agree that the scallop undergoing the shape-changes corresponding to the function 
Oi will not move on average, between the times t\ and t^. Likewise, the mollusk will 
be at the same place at the time £3 after performing either sequence corresponding 
to 0i or 62. One may also wonder where the animal would go asymptotically, as 
time goes to infinity. Following Purcell's reasoning, probably not very far and more 
precisely, exactly at the same distance as if the angle would have ranged from 1 to 
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Figure 2. Assume that the curves above give the value of the 
scallop's hinge angle 6 (see the left hand side of Fig. [I]) with respect 
to time. Although these shape-changes are neither time reversal 
invariant nor time periodic, the scallop is at the same position 
at the times t\ and t^. Besides, the displacement of the scallop 
between the times t = and t = £3 depends only on the values of 
0(0) and 0(ts) and not on the shape of the curve in between. 



over a finite time interval... because time does not matter in the low Reynolds 
world. This last property suggests that the hypotheses of the theorem could be 
restated in a purely geometric framework and one may even think at this point 
that, sticking to the scallop example, a reasonable statement could be something 
like: the displacement of the scallop is a continuous function of the angle range. As 
an obvious consequence, one would deduce that a bounded angle range implies a 
bounded displacement. This is true but unfortunately cannot be extended to the 
general case. Indeed, consider now an other example of swimmer, pictured on the 
right of Fig. [I] and called by Purcell the corkscrew (and whose way of swimming is 
quite obvious). The configuration space is the one dimensional torus S 1 and the 
rotation of the flagella is known to produce a net displacement of the hypothetic 
animal. On Fig. [3] is drawn the graph of a function giving the value of the angle of 
rotation, valued in S 1 , with respect to the time. It is not so easy to reiterate the 
exercise of the preceding example and to derive a purely geometric criterion (i.e. 
time independent) allowing to determine wether the displacement is bounded or 
not. The reason is that, by looking only at S 1 , it is not possible to determine how 
many tours have been performed by the flagella. To do so, we have to look at the 
angle as valued not in S 1 but in the universal cover of the manifold. The notion of 
universal cover will allow us to state a generalized and purely geometric version of 
the scallop theorem, which will be the second main contribution of the paper. 

1.3. Outline of the paper. In Section 2, we present an abstract framework and 
state a generalized scallop theorem for a class of shape-changing bodies, called linear 
swimmers. This is quite classical material, except for a topological interpretation 
of what a reciprocal motion is, which may be original. In Section 3, we prove that 
swimmers at low Reynolds numbers and high Reynolds numbers (in a potential 
fluid and with some symmetry assumptions) are linear swimmers and meet the 
requirements of our main theorem. Finally, a numerical simulation of a swimmer 
in a perfect fluid is given in Section 4. 
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Figure 3. The graph gives the angle of rotation of the 
corkscrew's flagella (pictured on the right hand side of Fig. [T]) val- 
ued in the one dimensional torus S 1 with respect to time. By 
looking only at the projection of the curve on S 1 , one cannot de- 
termine how many tours have been performed. This is an evidence 
that although Purcell's scallop theorem can be addressed using 
a purely geometric point of view, it cannot be done without the 
recourse to the notion of universal cover. 



2.1. General assumptions on the swimmer. We assume that any possible 
shape of the swimmer can be described by a so-called shape variable s living in 
a Banach space S (which can be infinite dimensional). So the shape-changes are 
described by means of a smooth shape function t G R+ \-> s(t) G S where t stands 
for the time and s = ds/dt is the rate of change. The variable q G Q, where Q is a 
smooth, finite dimensional Riemannian manifold, gives the position of the swimmer 
in the fluid. For instance, to describe the position of Purcell's scallop, we would 
choose Q = R because the scallop can only move along a straight line. Since its 
shape is thoroughly described by the angle we would have s = and S = R 
or S = R/27T as well. In [1 , we give an example of 2D-swimmer (see Fig. |4| in 
an infinite extent of perfect fluid with potential flow. In this case, the manifold Q 
is R 2 x SO (2), while the shape space S is an infinite dimensional Banach space, 
consisting of complex sequences s = (sk)k>i ( s k G C, fc > 1) and endowed with the 



Notice however that physical and mathematical constraints usually affect the pair 
(s, s) and lead to the definition of allowable shape function. It entails in particular 
that s is bound to remain in a subset of S and that s cannot take any value in 
S either. As an example, let us mention the constraint of self-propulsion, which 
means that although directly prescribed, the shape-changes have to result from 
the work of hypothetical internal forces, occurring within the swimmer (like for 
instance the work of muscles). This constraint prevents, for instance, translations 
to be considered as possible shape-changes. At this point, we also add the constraint 
that the path 7 = s(R+) be included in X, a one dimensional submanifold immersed 
in <S. From a physical point of view, it means that at any moment, there is only one 
degree of freedom in the shape-changes (see Fig. [5]) . We recall that two complete 
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Figure 4. Three examples of shapes of the authors' 2D- amoeba 
model described in [1 j. The frame is attached to the body and 
moves along with it. The shape variable s is here a complex se- 
quence and q = (r, R(0)) G R 2 x SO(2) where r is the position of 
the center of mass of the swimmer (expressed in a fixed Galilean 
frame) and R{0) is a rotation matrix of angle giving its orienta- 
tion. 




Figure 5. On the left, the path 7 is not included in a one di- 
mensional submanifold of S. It means that for both moments 
corresponding to s = si, there are two degrees of freedom for the 
shape-changes (each one corresponding to a branch). On the right, 
notice that the submanifold X can be bounded or not. 



examples of modeling (the low and the high Reynolds swimmers) are given in 
Section 03 

2.2. Dynamics of linear swimmers. Denote by T q Q the tangent space to Q at 
the point q, by TQ the tangent bundle to Q and by £(<S, TQ) the space of the linear 
mappings from S to TQ. We call linear swimmer, any model of shape-changing 
body whose dynamics has the form: 

(1) |q(t) = <F(q ) s),s(t)> ) (t > 0), 

where F : Q x S — >• £(<S,TQ) is a smooth function satisfying 

(i) (F(q, s), s) G T q Q for every q G Q and every s, s G S; 

(ii) There exists K > such that ||(F(q, s), s)|| Tq g < If ||s||s||s||s for every q G Q 
and every s, s G S. 
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Figure 6. To any path 7 in the configuration space S corresponds 
a path T in the space Q, once the initial point qo has been fixed. 
Because it is independent of the parameterization of 7, we can 
define the mapping s G 7 4 q(s) G T. Moreover, for any time- 
parameterization, the map s \-> q is linear. 



The Cauchy-Lipschitz theorem guarantees that, for any qo G Q and for any smooth 
allowable shape function s : R + — » S there exists a unique solution to ([!]) with 
Cauchy data q(0) = qo- 

2.3. Generalized scallop theorem. The main feature of linear swimmers' dy- 
namics is the following reparameterization property. 

Proposition 1. Let any allowable control function s : R + — >> S and any point 
qo G Q be given. Denote by q : R + — > Q the solution to Equation |Ip with initial 
condition qo- Then, for any C 1 function f3 : R + R +; the solution : R + — >> Q 
to Equation corresponding to the shape-changes := t G R+ i->> s(/3(t)) G <S ; 
initial condition q(/3(0)) satisfies q^ = q o /3. 

Proof The time derivatives of the functions qo/3 and q^ coincide, both being equal 
to (F(qy3, s^), &'q){}'. Since we also have qo /3(0) = q^(0), the conclusion follows 
from a direct application of the Cauchy-Lipschitz theorem. □ 

On Fig. [6] are presented some geometric interpretations of what a linear swimmer 
is. From Prop. [T] above, one can easily deduce: 

Proposition 2. For any T > 0, there exists a real number R > such that for 
any C 1 function /3 : R + —> [0, T] and for any initial condition qo G Q, the solution 
q/3 : R + — )> Q to Equation |Ip corresponding to the shape-changes t G R+ \-> 
s(/3(t)) G Q with initial condition qo ; remains in the ball of Q of center qo and 
radius R. 

Proof Fix qo G Q and denote by q the solution to Equation with Cauchy data 
qo- The interval [0, T] is compact and hence the set {q/3(R + ) : f3 G C X (R + , [0, T])} = 
{q(/3(R_i_)) : /3 G C 1 (R + , [0, T])} = q([0,T]) is also compact (because q is con- 
tinuous) and hence bounded by some constant, which in addition can be chosen 
independently of qo- Indeed, we have, for any t G [0, T]\ 

de(q(i),qo) < f ||(F(q( U ),s( M )),s( U ))|| Tq(u)Q d U < K f ||s( U )|| (S ||s( U )|| (S d M . 
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Since t G R+ ^ s(t) G S is smooth, the last integral is bounded for every t in [0, T] 
by J T ||s(u)|| 5 ||s(u)|| 5 dix < +00. □ 

Our topological version of Purcell's scallop theorem will be obtained by rein- 
terpreting Proposition |2j in the frame of differential geometry, using the classical 
notion of universal cover (see for instance [6] for an introduction to covering man- 
ifolds), which we now recall the definition: For any (finite dimensional) smooth 
connected Riemannian manifold X, the universal cover of X is a simply connected 
smooth Riemannian manifold X endowed with a canonical projection p : X — >> X 
enjoying the following property: For every x in X and y in X satisfying p(y) = 
there exists a neighborhood U y of y in X and a neighborhood U x of x in X such 
that p\u : U y —> U x be an isometric diffeomorphism. Any vector field v on X can 
be lifted to X by defining locally v(y) — (T y p)~ 1 v(p(y)) for any y in U y . Any curve 
r : [0,T] — )> X solution to the ODE f = v(r) can hence be lifted to X as well by 
choosing any base point yo in p _1 (r(0)), and by considering the solution to the 
ODE y = v(y) with initial condition y(0) = yo. 

The Banach structure of S induces a Riemannian structure gx on X. This Rie- 
mannian structure is compatible with the topology of X. Seen as a one dimensional 
manifold (endowed with its own topology), X can be either compact (or equiva- 
lently bounded for gx, and hence diffeomorphic to 5 1 , the one dimensional torus) 
or not (and hence diffeomorphic to R). In both cases, the universal cover of this 
manifold is R. With this material, we can restate Proposition [2] as follows: 

Theorem 1 (Generalized scallop theorem). Consider any smooth shape function 
t G R+ 1 — >► s(t) G X and any lift s : R+ — » X of s (this choice is unique up to the 
choice of the base point s(0) in p~ 1 (s(0))). If the subset s(R+) of X is of finite 
length, or equivalently if the topological closure in X o/s(R+) is compact, then any 
solution q : R + — ^ Q to Equation j7]) is bounded as well. 

On Fig. [7] and according to the theorem, only the shape-changes relating to the 
third case can result in locomotion. 

Proof. Assume that the path s(R+) C X is of finite length I (with I > 0) and 
denote by f : t G [0,1] \-> r{t) G s(R+) its arc- length parameterization. Then, there 
exists a smooth function (3 : R + — » [0,1] such that s = p o r o f3 and the conclusion 
follows from Proposition [2] □ 

The following comments are worth being considered: 

• The geometric hypothesis of the theorem is independent of the choice of 
the base point §0. 

• As already mentioned, the case where the shape function t G R+ \-> s(t) G S 
is not periodic and s(R+) of infinite length agrees with the hypothesis, 
whereas it is not covered by Purcell's original theorem. 

• The topological nature of Purcell's scallop theorem has been known for 
quite a long time. For instance, in [10], an interpretation of periodic shape- 
changes is given in term of retract. This result could be extended to non 
periodic and possibly non compact shape-changes by saying that the closure 
of s(R_|_) in X has to be homotopic to a compact set. 
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Figure 7. Denote 7 = s(R+). In cases 1 and 2, the shape-changes 
cannot result in locomotion. Locomotion is possible in the third 
case only. 



• In [2 j, the theorem is connected to the exactness of some closed differential 
1-form. Notice that in the simply connected universal cover, exactness and 
closedness of differential 1-form are actually equivalent. 

3. Swimmer at low and high Reynolds numbers 

In this Section we derive the Euler-Lagrange equations for low and high Reynolds 
swimmers. We show that, although the properties of the fluid are completely dif- 
ferent in both cases, the equations eventually agree with the general form of 
linear swimmers. In the modeling, we will assume that: 

(i) The swimmer is alone in the fluid and the fluid-swimmer system fills the whole 
space. It entails that all of the positions in the fluid are equivalent and the 
equations of motion can be written with respect to a frame attached to the 
swimmer. 

(ii) The buoyant force is neglected. 

(iii) The fluid-swimmer system is at rest at the initial time. 

3.1. Kinematics. The shape-changing body occupies a domain B of R 3 and T := 
R 3 \ 8 is the domain occupied by the surrounding fluid. We consider a Galilean 
fixed frame (ei,e 2 ,e 3 ) and a moving frame (e^e^e^) attached to the body. At 
any time there exists R G SO (3) such that e* = Rej and we assume that the 
origin of the latter frame coincides with the center of mass r G R 3 of the body. 
We introduce the notation q := (i?, r), which belongs to the Euclidean group Q := 
SO (3) x R 3 . The Eulerian rigid velocity field of the frame (e*) with respect to 
(e-j) is defined at any point x G R 3 by w r (x) := u) x [x — r) + v, where v := r 
and uj is the rotation vector defined by RR T x = uj x x for all x G R 3 . The shape 
changes are described by means of a set of diffeomorphisms % s , indexed by the 
shape variable s, and that map a reference domain (let say for instance the unit 
ball B) onto the domain B* of the body as seen by an observer attached to the 
moving frame (e*). The Eulerian velocity at any point x of the swimmer is the 
sum of the rigid velocity and the velocity of deformation: w = w r + where 
:= R(d s Xs(Xs(R T (% — r )) _1 )?s). It can be expressed in the moving frame: 
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Figure 8. Kinematics of the model: The Galilean frame (e^) and 
the moving frame (e*) with e* = Rej (R G SO (3)). Quantities 
are denoted with an asterisk when expressed in the moving frame. 
The domain of the body S* is the image of the unit ball B by 
a diffeomorphism Xs (s G S is the shape variable) and T is the 
domain of the fluid. The center of mass of the body is denoted r 
and v is its velocity. 

w* = w* + w*, where w* := u>* x (x*) T + v*, v* i? T v, x* R T (x - r), 
u)* : = R T u) and = (<9 s Xs(^*), s) (more generally, quantities will be denoted 
with an asterisk when expressed in the moving frame). The deformation tensor 
is F s := Dxs and, keeping the classical notation of Continuum Mechanics, we 
introduce J s := |det(F s )|. 

3.2. Dynamics. The density of the body can be deduced from a given constant 
density go > 0, defined in B, according to the conservation of mass principle: 
Q* ° Xs = Qo/ Js- The volume of the swimmer is Vol = J B J s (x)dx, its mass m = 
QoVol and its inertia tensor I* = g J B |Xs| 2 Id — Xs ^ Xsdx in (e*) and I s = M*i? T 
in (e-j). The deformations have to result from the work of internal forces within the 
body. It means that in the absence of fluid, the swimmer is not able to modify its 
linear and angular momenta. Assuming that the swimmer is at rest at some instant, 
we deduce that at any time J B (9 s Xs 7 s)dx = and f B (d s Xs,s) x Xsdx = 0. These 
equations have to be understood as constraints on the shape variable and will be 
termed subsequently the self-propulsion hypotheses. The fluid obeys, in the whole 
generality, to the Navier-Stokes equations for incompressible fluid: Qfjj^u— V-T = 
and V • u = in T for alH > (gf is the fluid's density, u the Eulerian velocity, 
D/Dt the convective derivative, T := 2/iD(u) -pld with D(u) := (l/2)(Vu + Vu T ) 
is the stress tensor and fi the dynamic viscosity). The rigid displacement of the body 
is governed by Newton's laws for linear and angular momenta: m^v = — J dB Tn da 
and ^(I 8 Lj) = — J dB Tn x (x — r) da (the rigid displacement is caused by the 
hydrodynamical forces only) where n is the unit vector to dB directed towards the 
interior of B. These equations have to be supplemented with boundary conditions 
on dB, which can be either u • n = w • n (slip boundary conditions) or u = w (no- 
slip boundary conditions) and with initial data: u(0) = uo, R(0) = Ro, r(0) = i*o, 
oj(0) — ujq and v(0) = vq. 
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We focus on two limit problems connecting to the value of the Reynolds number 
Re := gVL/ ji (V is the mean fluid velocity and L is a characteristic linear dimen- 
sion). The first case Re < 1 concerns low Reynolds swimmers like bacteria (or 
more generally so-called micro swimmers whose size is about ljim). For the second 
Re ^> 1, we will restrain our study to irrotational flows (i.e. V x u = 0) and so it 
is relevant for large animals swimming quite slowly, a case where vorticity can be 
neglected. 

3.3. Low Reynolds swimmers. For micro- swimmers, scientists agree that inertia 
(for both the fluid and the body) can be neglected in the dynamics. It means that 
in the modeling, we can set go = Qf — 0. In this case, the Navier-Stokes equations 
reduce to the steady Stokes equations —V • T = 0, V • u = and we choose no-slip 
boundary conditions u = w on dB. Introducing u*(x*) := R T u(Rx* + r) and 
p*(x*) = p(Rx + r), the equations keep the same form when expressed in the frame 
(e*, ej), namely: — V-T* = 0, V-u* = in J 7 * with boundary data: u* = w*. From 
a mathematical point of view, the main advantage is that the equations are now 
linear. Notice that since the equations are stationary, no initial data is required for 
the fluid. Newton's laws read f dB * T*nd<7 = and J dB * T*n x x*da = (it means 
that the system fluid- swimmer is in equilibrium at every moment. Indeed, since 
there is no mass, any force would produce an infinite acceleration). As already 
mentioned, the Stokes equations are linear. It entails that the solution (u*,p*) is 
linear with respect to the boundary data w* and we draw the same conclusion for 
the stress tensor a* because it is linear in (u*, p*). Observe now that w* is linear in 
the 3 components uj* (j = 1, 2, 3) of u?*, in the 3 components v* (j = 1,2,3) of v* 
and in s. We can then decompose any solution to the Stokes equations accordingly: 
u * = H 3 j= i UjUj +^* u i+3 + ( u 5> s), p* = Y? j= i u jPj+ v jPj+3 + (P*di s) and the stress 
tensor as well: T* = £? =1 u*. T* + v*T* +3 + (T*,s). Notice that the elementary 
solutions (u*,p*) as well as the elementary stress tensors T* depend on the shape 
variable s only. We next introduce the 6x6 matrix M r (s) whose entries M^-(s) are 
MT.(s) := J dB# e* • (T*n x x*)&cj = J dB# (x* x e*) -T*ndcr (1 < i < 3, 1 < j < 6) and 
M£(s) := f dB# e*_ 3 • T*nda (4 < i < 6, 1 < j < 6) and N(s), the linear continuous 
map from S into R 6 defined by (N(s), s) := {j QB# (T*, s)n x x*da, J QB# (T*, s)nda). 
We can rewrite Newton's laws as M r (s)q* + (N(s),s) = where q* := (u?*, v*) T G 
R 6 . Upon an integration by parts, we get the equivalent definition for he en- 
tries of the matrix MT(s): M[ j (s) := 2/i J B# D(u*) : D(u*)dx*, whence we de- 
duce that M r (s) is symmetric and positive definite. The same arguments for N(s) 
lead to the identity: ((N(s),s)) j = 2/i J B# D(u*) : L>((u*(s), s))dx*. We eventu- 
ally obtain the Euler-Lagrange equation: q* = — M r (s) _1 (N(s), s), or equivalently 
q = -R(q)M r (s)- 1 (N(s),s) where R(q) := diag(#,iT) and q := (u, v) T . Although 
this modeling is not new, the authors were not able to find the Euler-Lagrange 
equation in this particular form, allowing one in particular to deduce the following 
result: 

Proposition 3. The dynamics of a micro-swimmer is independent of the viscosity 
of the fluid. Or, in other words, the same shape changes produce the same rigid 
displacement, whatever the viscosity of the fluid is. 

Proof. Let (u*,p*) be an elementary solution (as defined in the modeling above) 
to the Stokes equations corresponding to a viscosity ji > 0, then (u*, (/i//i)p*) is 
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the same elementary solution corresponding to an other viscosity ft > 0. Since the 
Euler- Lagrange equation depends only on the Eulerian velocities u*, the proof is 
completed. □ 



3.4. High Reynolds swimmers. Assume now that the inertia is preponderant 
with respect to the viscous force (it is the case when Re <C 1). The Navier-Stokes 
equations simplify into the Euler equations: Qfj^u — V • T = 0, V • u = in J 7 
where T = — pld and we specify the boundary conditions to be: u • n = w • n on 
dB (slip boundary conditions). Like in the preceding Subsection, we will assume 
that at some instant, the fluid-body system is at rest. According to Kelvin's circu- 
lation theorem, if the flow is irrotational at some moment (i.e. V x u = 0) then, 
it has always been (and will always remain) irrotational. We can hence suppose 
that V x u = for all times and then, according to the Helmholtz decompo- 
sition, that there exists for all time t > a potential scalar function cp defined 
in J 7 , such that u = Vcp. The divergence- free condition leads to Atp = and 
the boundary condition reads: d n cp = w • n. Following our rule of notation, we 
introduce the function ip*(t,x*) := (p(t,R T (x — r)) (t > 0, G J 7 *), which is 
harmonic and satisfies d n ip* — w* • n on dB* . The potential ip* is linear in w*, 
so it can be decomposed into ip* = Y^j^^jPj + v j^-\-s + (^cb®) (^ ms P rocess 
is usually referred to as Kirchhoff's law). At this point, we do not invoke New- 
ton's laws to derive the Euler-Lagrange equation but rather use the formalism 
of Analytic Mechanics. Both approaches (Newton's laws of Classical Mechanics 
and the Least Action principle of Analytic Mechanics) are equivalent (as proved 
in [7]), but the latter is notably simpler and shorter. In the absence of buoyant 
force, the Lagrangian function C of the body-fluid system coincides with the ki- 
netic energy: C = m\\v*\ 2 + §u;* • + \ J B# g*\w* d \ 2 dx* + \ J T# Q f \u*\ 2 dx*. 
In this sum, one can identify, from the left to the right: the kinetic energy of 
the body connecting to the rigid motion (two first terms), the kinetic energy re- 
sulting from the deformations and the kinetic energy of the fluid. We can next 
compute that: J B * ^>*|w^| 2 dx* = J B qq\ (9 s x s , s) \ 2 dx (upon a change of variables) 
and fj-* ^>/|u*| 2 dx* = Qf\\7cp*\ 2 dx* . It leads us to introduce the so-called mass 
matrices Mj(s), whose entries (MJ)^-(s) are defined by (MJ)^-(s) := QfVtp* • 
V(p*dx* (1 < i,j < 6), and MT(s) := diag(I*, mid) +MJ(s). One easily checks that 
M r (s) is symmetric and positive definite. We define as well the linear map N(s) 
from S into R 6 by ((N(s),s)) j := QfVip* • V(^,s)dx* (1 < j < 6) and we can 
rewrite the kinetic energy of the fluid in the form: ^q* -Mj(s)q* + q* • (N(s), s). In- 
voking now the Least Action principle, we claim that the Euler-Lagrange equation 
is: SlC = where we have denoted Sl '= ~ ^ the Lagrangian differential 

operator connecting to the system of generalized coordinates (q, q). Introducing 
the impulses (n,P) T := M r (s)(u>*, v*) T and (A, L) T := (N(s),s) (homogeneous 
to momenta) and since (5z,q*,Q) = (O* x u>*,fT x v* — u>* x V*) T for any 
Q := (ft,V) T e R 6 (and Q* := (ft*,V*) T with ft* := R T fl, V* = i? T V), 
we deduce that (S L C,Q) = f t (U + A, P + L) • Q* + (n + A, P + L) • («* x 
u;*,f2* x v* - oj* x V*). The Euler-Lagrange equation is hence the system of 
ODEs: ^(n + A) = (n + A) xcj* + (P + L) x v* and ^(P + L) = (P + L) x w*. 
Since we have assumed that at some instant, the fluid-body system is at rest, we 
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deduce that (II + A) =0 and (P + L) = for all time (this solution is an ob- 
vious solution to the differential system), which can eventually be rewritten as: 
q* = -M r (s) _1 (N(s),s), or equivalently q = -R(q)M r (s)- 1 (N(s), s). 

3.5. Breaking the symmetry. Both models' responses to symmetry breaking 
are very different. If we assume that there is a rigid fixed obstacle in the fluid 
or that the fluid-swimmer system is confined in a bounded domain, the Euler- 
Lagrange equation for the low Reynolds swimmer still agrees with the form ([!]) 
and the scallop theorem still holds true. Indeed, although the matrix M r (s) is no 
longer independent of the position q and has to be rather denoted M r (s,q), the 
dynamics still reads q = — R(q)M r (s, q) _1 (N(s), s). Things begin to turn bad when 
additional degrees of freedom enter the game. It is the case when several swimmers 
are involved (this case is treated in [4]), when there is a moving rigid obstacle or 
when the swimmer is close to flexible walls. 

The high Reynolds swimmer is much more sensible to the relaxation of the 
hypotheses ( ipii ) and actually if any of these assumptions fails to be true, the 



Euler-Lagrange equation turns into a second order ODE containing a drift term. 
Obviously, the scallop theorem fails to apply in this case. We refer to [7] and [8] 
for details. 



4. When Purcell's scallop can swim... in a perfect fluid 

In Section [3j we have derived the Euler-Lagrange equations for both a low and 
high Reynolds swimmers. In both models, we have assumed that the fluid-body 
system was filling the whole space (the only boundary of the fluid was the one 
shared with the swimmer). In a potential flow, this hypothesis is necessary for the 
Euler-Lagrange equation to have the particular form required in the statement of 
Theorem [2] In this Section, we aim to show, through a numerical example, than 
the scallop theorem no longer holds true when the fluid contains in addition to the 
swimmer, a fixed obstacle. So, we consider the simple example of the scallop (as 
modeled in the original article of Purcell) swimming in a perfect fluid with irro- 
tational flow. Simulations have been realized with the Biohydrodynamics Matlab 
Toolbox (which is free, distributed under license GPL and can be downloaded at 
|http: //bht . gf orge . inria. f r/[ ). The scallop is made of two rigid ellipses linked 
together by a hinge. As shown on Fig[9| the animal is located close by a rectangular 
obstacle. This obstacle breaks the symmetry of the model as described in Section [3] 
and although the scallop is only able to flap, each stroke will make it get closer 



to the obstacle, until it eventually collides with it (see Fig. 10). A more physical 
explanation is that during a strokes, the fluid pressure (on which solely relies the 
hydrodynamical forces) is weaker between the obstacle and the scallop's left arm. 



5. Conclusion 

In this article, we have revisited Purcell's scallop theorem and proved that the 
common hypotheses on the sequence of shape-changes: time periodicity and time 
reversal invariance, although quite intuitive, are irrelevant from a mathematical 
point of view and have to be replaced by purely geometric conditions involving the 
universal cover of the configuration space. We have also shown that Purcell's result 
applies as well to swimmers at high Reynolds numbers and does not rely solely on 
the system's inertialess. 
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Figure 9. Purcell's scallop example (in a perfect fluid with poten- 
tial flow). The scallop is modeled as an articulated body consisting 
of two ellipses linked together by a hinge. The angle between the 
ellipses is a(t) = 7r/3cos(t) (t > 0). The flapping of the scallop 
does not produce locomotion (after completing a stroke, the scallop 
comes back to its exact initial position) in a fluid free of obstacles. 
However, this is no longer true in this example where a fixed im- 
mersed rigid solid breaks the symmetry of the model. The flapping 
motion generates a low pressure zone between the scallop's left arm 
and the obstacle, causing the animal to be attracted by the solid 
and eventually to collide with it. 




Figure 10. On the left hand side of the figure is plotted the in- 
coordinate of the center of mass of the scallop with respect to time 
and on the right hand side, the ^-coordinate of its velocity with 
respect to time. Owning to the presence of the fixed obstacle, the 
scallop undergoes a net displacement to the left. Notice that the 
velocity increases along with the number of strokes. 
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